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Abstract 

We propose a method to calculate the Greens function of a free massive scalar field 
on the lattice numerically to very high precision. For masses m < 2 (in lattice units) 
the massive Greens function can be expressed recursively in terms of the massless Greens 
function and just two additional mass-independent constants. 
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Lattice regularization offers a convenient tool to study non-perturbative features of field 
theories. The Lagrangian is formulated on a discrete space-time lattice and the corresponding 
finite path integral is evaluated numerically. In most instances, one is ultimately interested 
in the continuum limit by decreasing the lattice spacing to zero, but for very small lattice 
spacings the lattice path integral is plagued with ultraviolet divergences rendering numerical 
lattice simulations impossible. 

Perturbation theory plays an important role in providing the missing link between numerical 
simulations and the desired physical continuum limit. Perturbative calculations on the lattice 
are performed, e.g., to obtain renormalization constants of lattice operators or non-universal 
coefficients of /3-functions. If the physical situation contains separate scales, they are also 
necessary in order to match the perturbative short distance physics with the non-perturbative 
long distance physics. 

However, perturbative calculations on the lattice tend to be difficult, since the Feynman 
integrals are rather complicated functions of the involved momenta. In this respect, coordinate- 
space methods have been proven useful in the evaluation of Feynman diagrams and allow a very 
precise determination of the continuum limit of two- and even higher-loop integrals [1]. In [1] 
this technique was applied to massless propagators. The massive case was addressed in [3], 
but for space-time dimension less than four, and explicit results for the Greens functions were 
obtained only in one and two dimensions. In the present investigation, we extend the work of 
Liischer and Weisz [1] to massive scalar fields by presenting an efficient method to calculate the 
associated Greens function to very high precision. 

We will first generalize the recursion relation for the Greens function presented in [1] to 
massive propagators. The free propagator for a massive scalar field in position space reads in 
lattice units 




(1) 



where n — {rii, 712,113, n^) and 
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We define the forward and backward lattice derivatives 



V^G(n) = G{n + 11)- G{n) , V^G^ = G{n) - G{n - 11) 



(3) 



which yield the lattice Laplacian 
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with 
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The Greens function satisfies the equation 




(6) 
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from which one obtains after making use of hypercubic symmetry the identity 



8 + m2jG(0) -8G(^) = 1 (7) 
for any fi. The central identity utihzed in [1] is given by 

(v; + V,)GH=n^/7H , (8) 

where 

H{n) = r e'"-^ In f P + (9) 



is independent of /i. By ehminating H{n) via 



for p = X]^=i 7^ one obtains a recursion relation for the Greens function at different lattice 
sites ^ 

G{n + /i) = ^(n - /i) + ^ ^ ( [l + ^] G{n) - G{n - u)^ . (11) 

This recursion relation is at the heart of this method, as it allows to express the Greens function 
at any lattice site in terms of a few values close to the origin. Utilizing hypercubic symmetry 
of the lattice G{n) can be written as a linear combination of the values of the propagator at 
the corners of the unit hypercube 

^(0, 0, 0, 0) , ^(l, 0, 0, 0) , ^(l, 1, 0, 0) , ^(l, 1, 1, 0) , ^(l, 1, 1, 1). (12) 

The second value G(l, 0,0,0) can be eliminated by making use of the identity ((Tj). In the 
massless case, two more restrictions on the remaining four values of the Greens functions can 
be derived by reducing the relation in Eq. (jlll) to a one-dimensional recursion along the lattice 
axes and constructing two linear combinations of the Greens function values along the lattice 
axes which are independent of the lattice site n [1]. One can then perform the limit n — > cx) 
and compare it with the asymptotic form of the massless Greens function Gq yielding two more 
relations for the initial propagator values in (fT^ . Therefore, out of the five initial values three 
can be eliminated and the entire massless Greens function can be expressed in terms of just two 
values, say Go(0, 0, 0, 0) and Go{l, 1, 0, 0), which have been calculated by Liischer and Weisz to 
very high precision. 

For general non-vanishing mass, on the other hand, linear combinations of the propagator 
values along the axes which are independent of the lattice site cannot be constructed. (An 
exception to this are the imaginary masses = —4, —8, —12. For a similar observation in two 
dimensions cf. [3].) Therefore, with the method described above the initial values of the Greens 
function close to the origin cannot be reduced further by employing its asymptotic form. 

Nonetheless, there exists an additional constraint for the massive Greens function which 
allows us to express its values at any lattice site in terms of the massless Greens function 
Go and just two additional constants which are independent of the mass. As we will see. 
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these additional two coefficients can be calculated to very high precision yielding a precise 
determination of the Greens function for m < 2. This constraint is based on the identity for 
the modified Bessel functions /„ [4] 

A [/„_i(A) - J„+i(A)] = 2nh{\) . (13) 

The Greens function can be written as an integral over the Bessel functions 

G{n,,n2,ns,n^) = -j d\ e-"^'^/2-^"j„,(A)J„,(A)J„3(A)J„,(A) , (14) 

and utilizing the identity (fT!^ we obtain the relation 

G{ni + 1, n2, ns, rn) = G{ni - 1, ^2, ns, n^) + Go(ni + 1, n2, ns, ua) 

- Go{ni-l,n2,ns,n4) +ni dfi"^ G{ni,n2,n3,n^; fi^) , (15) 

Jo 

where fi is the mass of the Greens function in the integral and Gq is the massless Greens 
function. Therefore, knowledge of the functional dependence of G on the mass at lower lattice 
sites yields the mass-dependence at higher lattice sites. This constraint in combination with 
the recursion relation (|TT|l allows us to calculate G(l, 1, 0, 0), G{1, 1, 1, 0), G{1, 1, 1, 1) starting 
from the Greens function at the origin, the tadpole G{0, 0, 0, 0), if we make use of its functional 
dependence on m. As a matter of fact, it is well known that the tadpole can be expressed by 
the following small mass expansion [6,7] 

oo oo 

G(0, 0, 0,0) = ^ a, m^' + In ^ b, m^' (16) 

i=0 1=0 

with expansion coefficients Oj and 6j. We will show below that this expansion converges abso- 
lutely for masses < 4. Employing this expression, one uses then the relation 

(7(2,0,0,0) = G(0, 0,0,0) + Go(2, 0,0,0) -Go(0, 0,0,0) + / cZ/i^ G(l, 0, 0, 0; /i^) , (17) 

Jo 

where G(l, 0, 0, 0) can be directly derived from G(0, 0, 0, 0) via Eq. ((7j) and is easily integrated 
utilizing 

/ i/dix = — — — and / /iMn /i (i/i = — — — — (l + i)lnm^ — 1 . (18) 
Jo i + I Jq + 1) I J 

Next, the recursion relation (jllll can be employed to obtain G{1, 1, 0, 0) 

^(1, 1, 0, 0) = ^ ([8 + m^] G(l, 0, 0, 0) - G(2, 0, 0, 0) - G(0, 0, 0, 0)) . (19) 
Similar steps can be taken to obtain G(2, 1, 0, 0) 

G(2,1,0,0) = G(1, 0,0,0) + 6-0(2, 1,0,0) -6-0(1, 0,0,0)+ / rf/x^ G(l, 1, 0, 0; /x^) (20) 

Jo 
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which in turn leads to 



1 / r Tfi 1 \ 
1, 1, 0) = - ([4 + —J ^(1, 1, 0, 0) - ^(2, 1, 0, 0) - ^(1, 0, 0, 0)j . (21) 

Finally, one has 

G{2, 1, 1, 0) = 1, 0, 0) + 6-0(2, 1, 1, 0) - Go(l, 1, 0, 0) + / dfi^ G{1, 1, 1, 0; (22) 

Jo 

and 

G(l, 1, 1,1) = ^ ([8 + m^] G{1, 1, 1, 0) - 3 G(2, 1, 1, 0) - 3 G{1, 1, 0, 0)) . (23) 

We have thus calculated the Greens function at four lattice sites of the unit hypercube G{1, 0, 0, 0), 
G(l, 1, 0, 0), G(l, 1, 1, 0), G(l, 1, 1, 1) recursively by making use of Eqs. ©, (HH) and (jEI) which 
is sufficient to calculate G(ni, n2, na, 71,4) at all other lattice sites. As input we merely need the 
massless Greens function which is already known to high precision and the Greens function at 
the origin, i.e. the tadpole G(0, 0, 0, 0). 

For the calculation of the massive Greens function it is therefore sufficient to determine the 
coefficients Oj and bi in Eq. (jlfij) to very high precision. Moreover, one can express all expansion 
coefficients a, and bi just in terms of the leading coefficients 09,01,02 and 60 • Since oq is the 
massless tadpole Go(0, 0,0,0) and 60 can be calculated exactly, one is then left with the task 
to determine the constants 01,02 very precisely. To this end, consider the differential equation 
for the modified Bessel function /^(A) 

A2^/„(A) + A^/„(A) - (A^ + n') /„(A) = . (24) 
It follows that the product /o(A)^ of Bessel functions satisfies the differential equation 
A^:^ + lOA^-^ - 5X\AX^ - 5)-^ - 15A(8A2 - 1^ 



dX^ dX^ ' 'dX^ ' 'dX^ 

+(64A^ - 152A' + 1)-^ + 32A(4A' - l)Vo(A)' = . (25) 
dX J 

With the representation 

1 poo 

G(0, 0, 0, 0) = - / dX e-^"^/^-^^I,{Xt , (26) 
this translates into a differential equation for the tadpole [5] 

;8t^ - 160t^ + 512t) + (40t^ - 480t2 + 512) 



d{m?Y ' 'd{m?f 
+(50*» - 3040-^ + (15f^ - 32);^ + 0G(0. 0,0.0) = . (27) 

where t = 4 + m? /2. Inserting the ansatz (fTBj) into Eq. ()27|) each power m^* and m^Mnm^ 
must vanish separately. This leads to recursive relations between the expansion coefficients 



5 



tti, hi which are then expressible as hnear functions of the leading coefficients ao, Oi, 02 and 60 ■ 
The relations between the coefficients ^0,1,2,3 are given by 



3846o + 3072fei = , 

2O860 + 396861 + 1843262 = , 

626o + 251261 + 3110462 + 1105 9263 = , (28) 

while for the higher coefficients they read 

i%i-2 + (8 + 40i + + 80i^ + 40i^)6i_i + (832 + 2784i + 3632i^ + 2240i^ + 560i^)6i 
+ (15872 + 43008i + 43136i2 + 19200i^ + 3200^^)6^+1 + (73728 + 172032i 
+ 141312^2 + 49152i^ + 6144^^)6^+2 = (29) 

with i>2. The relations between the coefficients and 6, are 

2ao + 208ai + 3968a2 + 18432a3 + 6966o + 107526i + 43OO862 = 0, (30) 



Oo + 248ai + 10048a2 + 124416a3 + 442368a4 + 6OO60 + 190086i + 19968062 + 62668863 = 0, (31) 
and for i > 2 

i^ai_i + (8 + 40« + SOi^ + + 40i^)ai + (832 + 2784i + 3632^^ + 2240^^ + 560i^)ai+i 
+ (15872 + 43008i + 43136^^ + 19200i-'^ + 3200i^)ai+2 + (73728 + 172032i + 141312i^ 
+ 49152^3 + 6144i^)ai+3 + Ai^hi^2 + (40 + 160^ + 240^^ + 160i^)6i_i + (2784 + 7264i 
+ 6720^2 + 22AQi^)hi + (43008 + 86272? + 576002^ + 12800^^)6^+1 

+ (172032 + 282624i + 147456^^ + 24576i^)6i+2 = 0. (32) 

The constant oq has already been given in [2] with a precision of 10~^^® 

ao = 0.154933390231060214084837208107375088769161133645219 . . . , (33) 

and we are left with the precise determination of ai and 02- This can be accomplished by 
splitting the tadpole integral G(0, 0, 0, 0) into two parts [5] 

^(0, 0, 0, 0) = ^ / dX e-'^'^/^-^^IoiXy +^ dX e-'^'^'^-^^hiXf , (34) 
2 Jo ^ Ja 

where we leave the parameter A unspecified for the time being. The first integral is analytic 
in the mass and can be expanded in m^, while the logarithmic pieces proportional to Inm^ 
are contained in the second integral. The integration limit A is now chosen large enough, say 
A = 200 ~ 300, such that the integrand in the second integral can be replaced by the asymptotic 
form for the Bessel functions at large arguments. Hence we obtain 

1 /"A -1 ^ roo 

G{0, 0, 0, 0) = ^ / dX e-™'^/2-^^/o(A)^ + nJl'^n dX e-™'V2^-n (35) 
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with well-known expansion coefficients c„. The second integration can be performed analyt- 
ically, while the ffist integral can be evaluated numerically to very high precision, e.g., by 
employing a GauB-Legendre algorithm ^. Although the two single integrals depend on A, the 
dependence cancels out in the sum of both terms. Equation (jH^j) is therefore suited to determine 
the constants ao, cii, 02 and Bq. As a check we ffist calculate Oq 



ao 



Go(0, 0, 0, 0) = - / dX e-'%{Xr + -J^ . _ (36) 

"^0 n=2 ^ ^ 



and obtain agreement with the value given in Eq. ()33p up to a precision of 10~^^°. This precision 
is sufficient for our purposes, although it could be easily improved further by using slightly more 
computer time. The direct integration proposed here is thus an alternative method to the one 
advocated in [1] for obtaining the massless Greens function at the origin to very high precision. 
For the coefficients ai and 02 this procedure yields ^ 



n=2 

-0.030345755097111005404345333329331191871789937096381870985508127022 
928862366546240927128051477698201507815845611452036369186167716864948 
113458657909428..., (37) 



02 



n=2 

= 0.0027759274572839795869710267307723922239912302615036308120948317995 
820137296629367446462391602115470725514438235848939713163742942230032 
62761367623512 . . . , (38) 

where En is the exponential integral function defined by 

E4z) = f'-^dt, (39) 

and Tj denotes the i-th coefficient in the asymptotic expansion of En in m^. Utilizing the series 
representation of the exponential integral function 

En(z) = )~^^ [-\nz + ij(n)] - V ^'^^ , , , (40) 

m=0,m^n—l 

where ipif^) is digamma function defined by 

n-l 

m 



n—l _ 

^(l) = -7, ^(n) = -7 + V-, (41) 



m=l 



•^Since the integrand e " '*'*'/o(A)'* is an analytic function in A, the interpolation procedure converges. 
■^The first nine digits of oi can be found, e.g., in [8] 
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one can immediately read off the coefficients 



'AM A^ 



J — n 

n=3 



C3|*(3)-ln(^-j|-c,A-^c„— 



(42) 
(43) 



We have varied A within the range A = 200 . . . 400 and confirmed that the results are indeed 
independent of A up to a precision of 10"^^'^. The parameter Bq, on the other hand, originates 
solely from the second integral and can be calculated exactly: bo = l/167r^. From Eqs. 1)281 - 
it follows then that all coefficients bi can be calculated analytically. Having determined the 
constants Oq, ai and 02 to very high precision, we use the recursion relations for the coefficients 
and bi which follow from Eq. 1)77)) to calculate the last remaining value of the Greens function on 
the unit hypercube, G(0, 0, 0, 0). As an additional check, we have calculated the Greens function 
values G{0, 0, 0, 0), G(l, 1, 0, 0), G(l, 1, 1, 0), G(l, 1, 1, 1) by splitting the integral representation 
of Eq. (fT^ into two ranges [0, A] and [A, 00) analogous to Eq. and evaluating the integrals 
directly. The results obtained in this way agree with the ones achieved by applying the recursion 
relation up to 10^^'^^. With the method outlined above we thus confirm the form of the tadpole 
in Eq. 

Employing the recursion relations Eq. (f^H) and Eq. (jH^ we now illustrate that the small 
mass expansion of the tadpole (fT?)|) converges absolutely for < 4. For large enough i the 
recursion relation for the bi simplifies to 

bj+i _ _25 _ bi_3 _ 5 6^_2 _ 35 bi_i 
bi ~ 48 6144 6i 768 6^ 384 6, ' ^ ^ 

If we assume that for large i the ratio bi^i/bi converges asymptotically towards a constant value, 
then Eq. (jl^ has the solutions h^i/bi = ~j2^ ~Ts- have confirmed numerically 

that the ratio tends indeed towards — |. This implies that for large i the parameters bi scale 
as \bi\ oc 4^*. The logarithmic part in Eq. |T6|) converges thus absolutely for < 4. 

Assuming also a constant ratio aj+i/aj for large enough i and the boundedness of bi/ai we 
arrive at the same simplified recursion relation for the Oj 

di+i _ _25 _ CLi-s _ 5aj_2 _ 35aj_i , , 

~ ~48 ~ 6144 ~ 768 ai ~ 384 ^ ^ 

Again, we have checked numerically that both the ratio aj+i/aj tends towards and bi/ai is 
bounded. This implies that the polynomial part of the tadpole in fll6p and hence the entire 
small mass expansion of the tadpole converges absolutely for < 4. 

Next, we would like to confirm the results for the asymptotic behavior of the lattice Greens 
function for large \n\ = ^/nf + n| + n| + n|. In [9] Paladini and Sexton have derived an 
asymptotic form of G{n) for large \n\, Gas{n), in terms of modified Bessel functions of the second 
kind, Ki{m\n\). In the left panel of Fig. [T]we have plotted our results ln[G(n, 0, 0, 0)] for four 
different masses. We observe that with increasing n ln[G(n, 0, 0, 0)] becomes an asymptotically 
linear function in n. As for large n the Bessel functions Ki{mn) decrease exponentially as 
Kiimn) ~ \/2mn 6xp(— m n) , the different slopes in the plot reflect the different masses. In the 
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Figure 1: Left panel: Shown are the logarithms of the Greens functions G{n, 0, 0, 0) for masses 
= 0, ^, 1 (from top to bottom). Right panel: Ratio of Gas{n,0,0,0)/G{n, 0,0,0) for 

2 1 



right panel, the ratio Gasin, 0, 0, 0)/G{n, 0, 0, 0) is plotted for rn? = 1/2, demonstrating clearly 
the fast convergence of the asymptotic expansion towards the exact result. 

So far, we have discussed the case m < 2. The case m > 1, on the other hand, can be 
calculated in a straightforward manner. Starting from Eq. ()14|). we replace the Bessel functions 
by their Taylor series to obtain 



G{ni,n2,n3,ni) ^ 



1 r'^ °° 

- / dAe-'"'^/2-^"^rf,(ni,n2,n3,n4)A\ (46) 

-^0 i=0 



where the expansion coefficients di are derived from the product of the Taylor series of the 1^. 
Since the integrand in Eq. is absolutely convergent, summation and integration can be 
interchanged leading to 

G{ni,n2,n3,n4) = diini,n2, 713,714) h+ — j . (47) 

We have checked numerically that this series converges sufficiently fast for masses m > 1. E.g., 
for m = 1 we obtain a precision of 10~^^° by summing up the first n = 7000 terms in ()47|1 . while 
for = 8 the precision is 10~^^^ when taking the first 700 terms into account. Note, that 
the convergence behavior of the sum in Eq. ()47p worsens with decreasing m and is extremely 
inefficient for m <^ 1. It is therefore not suited to perform the limit m ^ 0. 

To conclude, we have presented the precise determination of the Greens function of a massive 
scalar field on the lattice which is based on a recursion relation for the values of the Greens 
function values at different lattice sites similar to the one obtained in the massless case [1]. 
We have shown that the Greens function for masses m < 2 is fixed by the knowledge of the 
massless Greens function and only two additional mass-independent constants Oi and 02, which 
can be determined to very high precision. We have also confirmed the usefulness of the large 
\n\ asymptotic expansion of G{n) given by Paladini and Sexton [9] and demonstrated its fast 
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convergence towards the exact result. For m > 1 the Greens function can be calculated directly 
by writing the Greens function as an integral over the Bessel functions and replacing these by 
the corresponding Taylor series. 

With minimal amount of computer time we achieve precisions of 10~^^° for the Greens 
function values at the origin. Such precise knowledge of the Greens function is expected to be 
sufficient for carrying out two- and higher-loop calculations. The applicability to higher-loop 
integrals is currently under investigation and could prove useful, e.g., in the recently proposed 
framework of lattice regularized chiral perturbation theory [10]. 

This work has been partially supported by the Deutsche Forschungsgemeinschaft. 
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